Vibrational spectroscopy of H^: 
hyperfine structure of two-photon transitions 



■ 



Jean-Philippe Karr, 1,2, Lj Franck Bielsa, 1 ' 2 Albane Douillet, 1 ' 2 Jofre 
Pedregosa Gutierrez, 1 ' 2 Vladimir I. Korobov, 3 and Laurent Hilico 1 ' 2 

1 Laboratoire Kastler Brossel, Universite Pierre et Marie Curie 

OO 

T12, Case 74, 4 place Jussieu, 75252 Paris, France 
2 Departement de Physique et Modelisation, Universite d'Evry Val d'Essonne 

S3 

£3 ' Boulevard F. Mitterrand, 91025 Evry, France 

■ 3 Joint Institute for Nuclear Research, 141980, Dubna, Russia 

(Dated: June 12, 2008) 

'Hh- Abstract 

We present the computation of two-photon transition spectra between ro-vibrational states of 
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I. INTRODUCTION 



The ion is the simplest stable molecule. It plays an important role, both as a bench- 
mark system for detailed studies of molecular energy levels jl| , and in astrophysics. However, 
there have been very few investigations concerning high- resolution spectroscopy of H^. Ra- 
diofrequency spectroscopy of the hyperfme structure has been performed on ions trapped 
in a Paul trap 2J. Rotational and ro- vibrational transitions close to the dissociation limit 
were investigated using microwave and laser spectroscopy on an ion beam 3j,|4|. The scarcity 
of experimental studies is mainly due to the fact that , being homonuclear, does not have 
a dipole-allowed rotational or vibrational spectrum (except in a small region close to the 
dissociation limit, where the lsa g and 2pa u electronic curves overlap). 

Two-photon ro-vibrational transitions are nevertheless allowed, and Doppler-free two- 
photon spectroscopy was proposed a few years ago as apromising new method for determi- 
nation of the electron-to-proton mass ratio m e /m p 0, 0|. Transition probabilities between 
L = states were computed in [7| (where L is the total orbital angular momentum quantum 
number), demonstrating the feasibility of two-photon spectroscopy using Av = 1 transitions 
around 8-12 /jm. Among these, transitions lying in the spectral range of CO2 lasers (9-10 
fim) are especially attractive for frequency metrology, because of their high output power 
and stability. Even if there is no coincidence of Ht transitions with the CO? lines, a CO2 

n 

laser can be used as a frequency reference for a tunable quantum cascade laser (QCL) [8j. 
We have built an experiment designed to probe the (v — 0, L — 2) — > (v' = 1,L' = 2) 
transition at 9.166 jj,m 

The aim of this paper is to present the computation of two-photon transition spectra in 
H2", including the effects of hyperfme structure. In Sec. [Til the hyperfme effective Hamil- 
tonian obtained in jlOj is diagonalized, and the hyperfme states are written explicitly. In 
Sec. PH we recall the theory of two-photon transition probabilities. The transition matrix 
elements between hyperfme states are expressed as a function of reduced matrix elements 
involving only the orbital part of the wave functions, which are calculated using the same 



variational method as in 



101 ]. In order to avoid huge data, only the spectra of the transi- 



tions (v = 0,L) — > (v' = l,L) with < L < 3 are presented 11|. One reason for this choice 
is that the hyperfme structure is essentially determined by the value of L. Moreover, 
the considered L values are the only one which are significantly populated when ions 
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are created by electron impact ionization on H2 at room temperature, and the frequencies 
of these transitions are sufficiently close to a CO2 line to allow their excitation by a laser 
system discussed below. 



II. HYPERFINE STRUCTURE OF 
A. Hyperfine Hamiltonian 

The following notations are used throughout this paper: S e and Ii, I2 are respectively 
the electron spin and the spins of both protons, with S e = 1% = I2 = 1/2. We introduce 
the total nuclear spin I = Ii + 12, where I is equal to or 1. The total orbital angular 
momentum quantum number is denoted L. Note that due to the Pauli symmetrization, and 
taking into account that the electron is in the ground lsa g state, the total nuclear spin I is 
equal to when L is even, and to 1 when L is odd. 



The hyperfine effective hamiltonian of the molecular ion is taken in a form 101 ]: 



H Ms = Ml-S e )+c e (L-S e )+c 7 (L-I) 
di f2 r2 . 



+ (2L _ lK2L + 3 ) U L (I • Se) - [(L • I)(L ■ Se) + (L • Se)(L • I} 

d 2 / 1 T 2 T 2 1 



+ 7 77 7 -L I — -(L • I) — (L • I) . (1) 

(2L - l)(2L + 3) V3 2 y ; y J J K ' 

The numerical values of the coefficients bp, c e , cj, di, d,2 have been computed with a relative 
accuracy of 0(a 2 ) [10] using a variational method, for all ro-vibrational levels (v, L) with 
< L < 4 and < v < 4. 

If / 7^ 0, the strongest coupling is the spin-spin electron-proton interaction, i.e. the first 
term in equation ([T]). This interaction determines the principal splitting of the ro-vibrational 
levels of H2". With this consideration in mind, the preferable coupling scheme of angular 
momentum operators is 

F = S e + I, J = L + F. (2) 

The possible values of F and J, as well as the number of hyperfine levels, are given in 
Tabled] for each value of L. The hyperfine structure is much simpler for the states of even 
L, where only the value F = 1/2 is allowed since the total nuclear spin is zero. 
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B. Hyperfine states 



In order to obtain the hyperfine eigenstates and frequency shifts, it is necessary to diag- 
onalize the Hamiltonian ([T]). This is immediate when L is even: the effective Hamiltonian 
reduces to c e (L • S e ), and can be written 



hfs 



2 



(J 2 



(3) 



Its eigenstates are the states \v,L,S e = |,I = 0, F — |, J, Mj) coupled according to the 
angular summation scheme (j2J), the corresponding energy shifts are : 

L + l 



1 1 T 1 
v, L, -, 0, -, L — 
''22 2 

T 1 1 T 1 
v,L,—, 0, -, L + — 
''22 2 



-f^hfs v, L, -, 0, -, L-- 

1 1 1 

Fhf s |?;,L, -,0, -,L+- 



L 



(L ? 0), (4) 
(5) 



All energy shifts for L = 0, 2 and u = 0, 1 are given in Table [TT1 The relative theoretical 
accuracy is 0(a 2 ), corresponding to the limit of the Breit-Pauli Hamiltonian used in 10| to 
compute the hyperfine coefficients. The numerical accuracy is higher, which is why more 
digits are given here, as well as in Table II I II below; the extra digits will become useful when 
higher-order corrections to the hyperfine structure are computed. 



The case of odd L is more complicated. The operators involved in the expression of H^s 
are I • S e , L • S e , L • I, L 2 and I 2 . Note that they all commute with L 2 , S 2 , I 2 , J 2 and J z , but 
the terms L • S e and L • I do not commute with F 2 . As a consequence, F is an approximate 
quantum number only. There is a degeneracy in Mj, so that it suffices to diagonalize the 
restriction of to a subspace of given Mj. In the following, Mj is set to 1/2. 
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Se 
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n 
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even 


1 
2 





1 

2 
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odd 


1 
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1 

2 

3 
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L— |, L— 2 , £+§ 


6 



TABLE I: Possible values of F and J as a function of L. n is the number of hyperfine levels. 
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L 


V 


J = L-l/2 


J = L+l/2 










0.0000 





1 




0.0000 


2 





-63.2438 


42.1625 


2 


1 


-59.3574 


39.5716 



TABLE II: Hyperfine splitting (in MHz) for the ro-vibrational levels (v,L) with L = 0, 2 and 
v = 0, 1. All digits are converged. The relative theoretical accuracy is 0(a 2 ), which corresponds 
to an uncertainty of a few kHz. 



Let us consider a set of states: 



\ F = U = L + I), 



J = L+\), \F = l,J=L+l), 



2i" 1 2" 1 2' 2" 1 2' 1 2" 

|-F = §, J = L — |), |F=|,J = L — |), |F=|,J = L— |), 

where the last ket exists only if L > 3. We will refer to them as pure states. The matrix 
representing H h{s in this basis can be derived by use of the following relations 



I S e 

(FJ\L-S e \F'J) 

(FJ\L-I\F'J) 
L 2 I 2 



1 



f 



\J+L+F 



1 

" 2 

L 1 L 
F' J F 

L 1 L 



11 

T 



(6) 



L(L + 1)(2L + 1) (S e , I, F\\S e \\S e , I, F') (7) 



-l) J + L + F l ~ ' }\/L(L+l)(2L+l)(S e ,I,F\\I\\S e ,I,F') (8) 
F' J F 



2L(L + 1) 



(9) 



and the reduced matrices of S e and I on the subspaces S = |F = | , F = || (see Eq. (91) 
of Ref. Q) 



>/T5 


2 


3 


V3 


2 


V6 




6 



/ 2y / 15 2 

3 V3 

2 2y^6 

V 3 



(10) 
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Since there is no coupling between different J states, the shape of H^s is the following: 



Ms 



^ 

B C 

C D 

E G 

G H 

K 



The nonzero coefficients are calculated from equations (!6l fTU|) : 

1 2dt + d 2 \ 



A 
B 

C 
D 
E 

G 
H 
K 



1+ 



L 
2 
L 



c e + 2cr — 
-3 



6 



(c e + 2c/ 



L(2L + 3) 



(c e - Cj) - 



3 2L + 3 / 

L + 3 2di + d 2 
~6 2L + 3 



L . 

-OF - 7T [C e - 4C/J 

o 

y - -g- (ce + 2 CJ ) + 



6V2L + 3 
L-22d l + d 2 



(di - d 2 ) 



6 2L 



(L + 1)(2L- 1) 



-Of H ; — [C, 



{C e - Ci) 

4c z ) 



6V2L - 1 



(rfi - da) 



6 

6 F L + l ( , _ , 12d! + d 2 ^ 

c e + 2c/ H 

2 2 I 3 2L - 1 



(11) 



(12) 
(13) 

(14) 

(15) 
(16) 

(17) 
(18) 
(19) 



The eigenstates of J = L ± | are pure states of angular coupling: |d , L, S e = 1/2, / = 1, F — 
3/2, J = L ± 3/2), while the eigenstates of J = L ± | are linear combinations of F = l/2 and 
F = 3/2 states, obtained by diagonalization of the 2x2 sub-matrices appearing in ( fTTT) : 

|«,L > 5 e ,J,F,J=L±i)=Cf|i;,L,i,l,i > L±|) + | u, L, |, 1, |, L ± |). (20) 

We will refer to them as mixed states. The coefficients Cf and are calculated in Table ITO1 
together with the hyperfine frequency shifts, for L = 1, 3 and v — 0, 1. The mixing between 
F = 1/2 and F = 3/2 states is weak, so that the states can be labeled by the dominant F, 
noted as F. The hyperfine splitting of the first ro- vibrational levels is shown in Figs. [1] 
and [3 



6 



v=l 



v=0 



■J=l/2 





L=0 



J=5/2 

105.4 MHz 
J=3/2 



L=2 



FIG. 1: Hyperfine splitting of the ro-vibrational levels (v,L) with L = 0,2 and v = 0,1. The 
spacings between hyperfine states are proportional to the frequency difference. That scale is not 
respected for the rotational and vibrational spacings. 



v=l 



v=0 



449.4 
" F=3~/2 " 



J=3/2 
1=512 



-o»o.o J=1/2 

F=l/2 J=3/2 

461.5 ^^^^ J=3/2 
F=3~/2 ~ ' J=5/2 
J=l/2 

-923.0 == J=I/2 

F=l/2 J=3/2 

I. ! 





L=3 



FIG. 2: Same as Fig. [H with L = 1, 3. Within a given i 7 multiplet, the spacings between J states 
are proportional to the frequency difference. That scale is not respected for the other spacings. 
The frequency shift of the center of the F multiplets with respect to the spin-independent level, 
are indicated in MHz. 
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III. TWO-PHOTON TRANSITIONS 



A. Two-photon transition operator 

In this paragraph, we present the general theory of two-photon transitions with arbitrary 



excitation polarizations, as developed by G. Grynberg in [13|. Let us consider an Hj^ ion 
irradiated by two beams of polarizations E\ and e 2 . The transition probability between two 
states |0) and by absorption of one photon in each wave is proportional to 

M S Q^M\ 2 (21) 

where 

Qe±e 2 = 2 (Qtl£2 "I" Qc2Cl) (22) 

is the two-photon transition operator, with 

Q ei£2 = d- ei -^l_d-e 2 (23) 

In this expression, d is the dipole operator, H is the full hamiltonian and E the intermediate 
state energy. If the excitation polarizations are chosen among the standard polarizations tt, 
a + and <x~, the two-photon transition operator reads 

Qqiq2 = 2 (Qqi92 Qq2Qi) ) Qqiq2 = ^qi _ < ^ < ? 2 

where d qi (<& = —1,0, 1) are the standard components of d. Tensor Q qiq2 has a rank 2 and 
can be represented in terms of irreducible tensors: 



Q[ k) = £<fc?|H<7i?2> fc=0,l,2. (25) 

Inverting this expression, one finds: 

S Qq,qi= £ a?Q?+$ ] Q$\ (26) 

q=-2 

where 

= (ll gi g 2 |A:g). (27) 

Table IIVI gives the values of the coefficients a q ® for all combinations of the standard polar- 
izations. 
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B. Two-photon matrix elements between hyperfine levels 



We consider a two-photon transition between the hyperfine states \<f>) = 
\v, L, S e , I, F, J, Mj) = \g, J, Mj) and \t/j) = \v' , L', S e , F , F', J', M'J) = |e, J', M'j) with stan- 
dard excitation polarizations qi, q^. In order to simplify the expressions, we restrict the 
presentation to a case where the initial and final states are pure states; the results will be 
generalized at the end of the paragraph. Using the Wigner-Eckart theorem, the two-photon 
matrix element between \<f>) and \ip) may be expressed as: 

[gJ\\QW\\eJ>) 



Q**W) = E<° (J'kM' jq \JMj) ^ L q = qi + q2 . (28) 



A- 



The states \g, J,Mj) and |e, J 1 , M'j) are degenerate in Mj or M'j. If we assume the initial 
level to be unpolarized, the transition probability averaged over Mj and M'j is proportional 
to the averaged squared matrix element 



S Q qi A9J^eJ')] 2 = -^— £ \(g,J,Mj\ s Q qiq2 \e } J',M:j)\ 2 . (29) 

ZJ 1 Mj,M'j 



n 

Using the orthogonality relations of Clebsch-Gordan coefficients [12J, one obtains [1 



2 



,2 1 (gJ\\Q {k) \\eJ') 
The dipole operator and hence acts on the spatial variables only. Using Eq. (89) of 



Ref. 



12j one can write the reduced matrix elements of using only orbital wave functions: 



gJ\\Q {k) \\eJ') = V^,fK-1) J ' + ^ F ^v / 2J+T v / 2J^| ^ * ^ j (vL\\Q^\\v'L') . (31) 

In the case where the initial and final states are not pure basis states (i.e. for odd L and 
J = L ± 1/2), they can be written according to Eq. (|20|) : 

\g,J)= £ C Fi \v,L, 1/2, l,Fi, J) (32) 

P. -I 3 
2'2 

This expression can also be applied to pure states, where one coefficient is equal to zero 
and the other is equal to one. It is then straightforward to generalize Eq. ( 13T1) : 

gJ\\Q^\ef) = Sjj, Y, 6 FuF ,C Fi ^ L k L \ (vL\\Q^\\v'L' 

' F iiF > 33 [J'Fi J J X 

(33) 



C. Selection rules 



Since the two-photon transition operator is a sum of operators of rank and 2, the 
states |0) = \v, L, S e , I, F, J, Mj) and \ip) = \v', L', S e , I', F', J', M'j) can be coupled only if 
\L -L'\<2 and \ J - J'\ < 2. 

For the ro- vibrational states of Hj, the total nuclear spin is / = when L is even and 
1=1 when L is odd. The two-photon transition operator acts on the orbital variables only, 
which explains the 5n' factor in Eq. (I3~TI) and gives the selection rule AL = 0, or ±2. 

For the same reason, if we consider pure states, we get the selection rule AF = 0, as can 
be seen from the Sff' factor in Eq. (1311) . However, due to the mixing between F = l/2 and 
F = 3/2, transitions between mixed states of different F are weakly allowed. 

There are additional selection rules on Mj and M'j depending on the beam polarizations 
as can be seen in Table HVl A difference M'j — Mj must be equal to for tctt and cr + cr_, +1 
(— 1) for 7ta + (7rcr_) and +2 (—2) for a + a + (er_er_). 



D. Reduced orbital two-photon matrix elements 

The last step consists in the numerical computation of the reduced matrix elements 
vL\\Q^\\v'L'\. This is achieved using the variational approach outlined in Ref. [h]]. Briefly, 
the wave function for a state with a total orbital angular momentum L and of a total spatial 
parity tt = (— 1) L is expanded as follows: 

(34) 

G% 2 (R,n,r 2 ) = J2 {C n Re[e- a " R -^ ri -^} + D n lm[e- a " R -^ ri -^}}. 

n=X 

where the complex exponents a, (3, 7, are generated in a pseudorandom way. The use 
of complex exponents instead of real ones allows to reproduce the oscillatory behavior of 
the vibrational part of the wave function and improves the convergence rate. Since very 
high accuracy is not required for transition probabilities, relatively small basis lengths of 
iV = 700 — 1000 were used, providing a relative accuracy of a few parts in 10 9 for the 
nonrelativistic energies, and a few parts in 10 5 for the matrix elements. 

The reduced matrix elements (^vL\\Q^ \\v'L'\ are divided into three terms corresponding 
to the possible values L — 1, L+l, L for the angular momentum of the intermediate state. 
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The three following terms are evaluated numerically: 

(vL\\d\\v"L-l) (v"L-l\\d\\v'L') 



a 



-E 

v" 

-E 

v" 

E 



{2L+1){2L' + 1){lo- E^,^) 
(vL\\d\\v"L+l) <^ // Zv + 1| l-^'Z.') 



(2L+1)(2L' + 1)(lu- E V , /LH ) 
(vL\\d\\v"L) (v"L\\d\\v'L') 



(35) 
(36) 
(37) 



7' ^{2L+l){2L'+l){u - E V ,, L ) 

where E v ii^ h n is the energy of the intermediate state \v"L") and uj = [E v il' — E v l)/2 is the 
photon energy. The reduced matrix elements of Q^ k > are related to a_, a+, a in the following 
way: 



vL 



v'L 



V2L + 1 
vL\\QW\\v'L-2 



a + a+; 



V2L + 1 
uL||Q( 2 )||i;'L 



'2L-3 
2L-1 



V2T+1 
i;L||Q( 2 )||t/L + 2 
+ 1 



v/6 



(2L+3)(2L-1)L(L+1) 



a 



+ 



L(2L-1) L(L+1) (2L + 3)(L+1 



'2L+5 
2L+3 



The reduced matrix elements of and for the transitions (v = 0,L) — > (t>' = 1,L) 
with < L < 3, are given in Table PVT. 



E. Two-photon transition spectra 

Using Eqs. fl30l) . (1331 . the orbital reduced matrix elements given in Table [V] and the 
mixing coefficients given in Table 1111} we have computed the two-photon matrix elements 
for the four transitions (v — 0, L) — > (v' — l, L) with < L < 3. They are given in Tables IVTl 
IVIII and IVIIII for L = and 2, L — 1 and L = 3 respectively. The corresponding spectra, for 
three different choices of standard polarizations: linear-linear, a + a + and cr + cr~, are shown 
in Figs. El land El 

Simplest is, of course, the L = case, where is no hyperfine splitting. The transition 
probability had been computed in Ref. |7| for linear-linear polarizations. Note that the 
transition is forbidden for cr + <j + polarizations, because of the selection rule AMj = 2. In the 
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20 



-i 1 1 1 1 >f 1 1— 

-60 -40 -20 20 

v-v (MHz) 



40 



40 



60 



60 



S Q 



91,92 



in atomic units between the ro- 



FIG. 3: Averaged two-photon matrix elements A 
vibrational levels {v = 0,L) and (v = 1,L) with L = 0, 2 (from Table IVT|) . The spectrum is 
centered around the spin- independent transition frequency given in Table IIXI 

L = 2 case, there are two intense AJ = lines shifted by a few MHz and two weak AJ = ±1 
lines shifted by about 50 MHz. 

The spectra are more complex for odd values of L. They consist in one main cluster of 
intense AF = lines which is about 50-100 MHz wide, and two satellite clusters of very weak 
lines (corresponding to AF = ±1) about 600-700 MHz away. The total number of lines is 
25(34) for L — 1(3) but the most intense are those of AF = AJ = 0; there are 5(6) of them 
for L = 1(3). 

Whatever the value of L, all the favored transitions are between states with similar 
spin structure (i.e. same values of F,J). This feature makes them especially attractive 
for metrological purposes. Indeed, in such pairs of homologous hyperfme states, systematic 
shifts like the Zeeman shift (see Ref. [ljjl) are expected to have similar values, so that 
the shift of the transition frequency will be much smaller. The same is true for hyperfine 
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FIG. 4: Same as Fig. with L = 1. 

structure corrections to the transition frequency: as can be seen e.g. from Figs. H] and 0, 
the most intense lines span a frequency interval of less than 25 MHz because the spin- 
dependent corrections to the initial and final state energies partially cancel each other. For 
this reason, the theoretical uncertainty on the frequency of these transitions is much smaller 
with respect to the other ones. On the whole, the favored transitions benefit at the same 
time from a smaller sensitivity to systematic effects, and from potentially more accurate 
theoretical predictions. 



F. Orders of magnitude 



The two-photon transition probability at resonance is 



'47rof\ 2 4 
he I T f 



<?1,<?2 



(42) 
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FIG. 5: Same as Fig. with L = 3. 

where a is the Bohr radius, Tf the instrumental width of the transition, and / is the laser 
beam intensity. The above results show that the averaged two-photon matrix element for the 
favored transitions does not depend critically on the value of L, F, J but strongly depends 

r„ - 1 2 

on the excitation polarizations; we have the typical values Q qi) q 2 ~ 0.2 for the case of 

r„ - -|2 

linear-linear or cr + a_ polarizations and Q qi , q2 ~ 0.02 for the <j+<j+ case. In this subsection 
we evaluate the two-photon transition probability, using the parameters of our experiment. 
Our excitation source is a QCL phase-.oched to a C0 2 .a Ser fl, which delivers a hnea riy 
polarized beam with a cw power of about 90 mW. A Fabry Perot cavity of finesse 1000 is 
built around the ion cloud. The QCL requires a strong optical isolation due to its extreme 
sensitivity to optical feedback from the high finesse Fabry Perot cavity. An optical isolation 
of more than 23dB (with 90% transmission) can be achieved using an optical diode made of 
a grid polariser and a quarter-wave plate, which implies working with 0" + cr + polarizations; 
the isolation ratio is limited by the polarizer extinction ratio [15] . An additional isolation of 
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A.O.M. 

n 



FIG. 6: Setup of our experiment for excitation with circular polarizations. A.O.M. stands for 
acousto-optic modulator, P for polarizer, CD for cold detector. 

6 dB is obtained using an acousto-optic modulator with a polarization-dependent efficiency. 
The setup we have implemented is shown in Figure El The overall transmission of those 
optical elements including the alignment mirrors is 60% so that 54 mW of optical power are 
injected into the high finesse cavity. The transmitted power at resonance is about 10 mW; 
from transmission and reflectivity measurements we estimate the mirror transmission and 
losses to about 0.001, so that the incident power on the ions is then P ~ 10W in a 
beam of waist Wq = 1mm. The intensity on the beam axis is IPli^wi ~ 6.4 W.mm" 2 . The 
instrumental width is essentially the laser width Tf = 2irx 2.6 kHz [16j, since the width of 
the excited state is extremely small, all the ro-vibrational states of being metastable. 
From equation (l4"2l) one obtains a transition probability T ~ 0.7 s" 1 with <y + a + polarizations. 

The transition probabilities are higher by about one order of magnitude in the linear- 
linear polarization case. Due to this, even with a circularly polarized beam it is still more 
advantageous to probe the AMj = transitions in a transverse magnetic field, which must 
be sufficiently strong to separate the three components AMj = 0, ±2. A field in the 100 
mG - 1 G range is enough, as estimated in [ijj]. The incident intensity is then decomposed 
into 50%, 25% and 25% of linear, <r_ and cr + polarizations respectively. A factor of 2 is 
lost on the excitation beam intensity (hence 4 on the transition probability), but this is 
more than compensated by the difference in the two-photon matrix element. With these 
parameters, the transition probability is T ~ 1.7 s -1 , a large enough value to observe a 
two-photon transition in Paul traps where the ion lifetime is typically of several seconds. 
Further improvement can be achieved either by a tighter focusing of the laser (and a smaller 
ion cloud section in order to minimize transit-time broadening) or by reducing the laser 
linewidth. 

A factor of 4 on the transition probability can be gained by using a linearly polarized 
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p k 
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A/4 A.O.M 



CD 



PZT 




FIG. 7: Proposed experimental setup for excitation with linear polarizations. A.O.M. stands for 
acousto-optic modulator, P for polarizer, CD for cold detector. 

excitation beam. In this case, the standard optical isolation technique relies on Faraday 
isolators. In the 9 micron range, a 45° polarization rotation with reasonable magnetic 



fields can only be o 
insertion losses 17 



jtained using n-doped InSb wafers under cryogenic conditions, with high 
. In addition, to our knowledge they are no longer commercially 



available. Other ways of achieving isolation must be sought. Figure [7] shows our proposal 
for a high-transmission high-isolation device for linear polarization. It takes advantage of 
both constructive and destructive interference by a Fabry Perot cavity. 

A Fabry Perot cavity of free spectral range 4/ is locked on resonance with the laser of 
frequency i>l, resulting in a high transmission. The transmitted beam is frequency shifted 
to vl + / by an acousto-optic modulator driven at a frequency /, and injected into the high 
finesse Fabry Perot cavity surrounding the ion cloud. On the way back to the QCL, the 
reflected beam is diffracted again by the acousto-optic modulator and shifted to u L + 2f . It is 
then exactly off resonance with the first Fabry Perot cavity that provides optical isolation. 
To summarize, the high transmission is due to constructive interference and isolation to 
destructive interference. Optical isolation of the QCL against the feedback from the first 
cavity can be achieved using an optical diode as discussed above. A second quarter-wave 
plate turns the polarization back to linear at the output of the isolation cavity. 

The performances of this setup can be estimated as follows. The transmission at resonance 
for a Fabry Perot cavity made of two identical mirrors of reflectivity R, transmission T and 
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losses P with R + T + P = 1, is 



Tc - [l+P/(l-i2-P)] 2 ' (43) 



and the off resonance isolation ratio expressed in dB is given by 

;i -R-Pf 



-10 log 



10 



(44) 



(1 + R) 2 

Using low-losses mirrors with R =0.98 and P =0.001, one obtains T cav =0.9 and I =40 dB. 

IV. CONCLUSION 

We have presented a derivation of the hyperfine structure of two-photon transition spectra 
in the molecular ion, and applied it to several rotational components of the fundamental 
vibrational transition (v — 0,L) —>■ (v' = l,L). It was shown that the most intense lines are 
those between pairs of homologous hyperfine states (v, L, F, J) — > (V, L, F, J). Our estimate 
reveals that observation of such lines in Doppler-free spectroscopy is feasible with present- 
day laser sources. We have also proposed an experimental setup allowing to probe the two- 
photon transitions with linear-linear polarizations. Let us point out that the experimental 
task of find™ the transition frequency is made easier by recent progress in theoretical 
predictions [19j. The current theoretical uncertainty on the spin- independent frequencies 
(given in Table IPX)) is about 13 kHz [^(J, while the uncertainty due to hyperfine corrections 
is of the order of 5 kHz, due to partial cancellation between the shifts of initial and final 
states. We have also shown that such transitions have a very low sensitivity to external 
magnetic fields jl4]. 
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L 


V 


F 


J 


AE'hfs 




ct,cf 








3/2 


5/2 


474.1063 




[0,1] 








3/2 


3/2 


481.9534 


[0.015612,0.999878] 


1 





1/2 


3/2 


-930.4332 


[-0.999878,0.015612] 






3/2 


1/2 


385.3985 


[0.038891,0.999243] 






1/2 


1/2 


-910.7579 


[-0.999243,0.038891] 






3/2 


5/2 


461.2574 




[0,1] 








3/2 


3/2 


468.5247 


[0.015074,0.99 


9886] 


1 


1 


1/2 


3/2 


-905.7836 


[-0.999886, 0.015074] 






3/2 


1/2 


377.9948 


[0.037345, 0.999302] 






1/2 


1/2 


-887.2491 


[-0.999302, 0.037345] 






3/2 


9/2 


507.2568 




[0,1] 








3/2 


7/2 


489.5257 


[0.042115,0.999113] 


3 





1/2 


7/2 


-941.1034 


[-0.999113,0.042115] 






3/2 


5/2 


423.6342 


[0.061812,0.998088] 






1/2 


5/2 


-894.6614 


[-0.998088,0.061812] 






3/2 


3/2 


341.5540 




[0,1] 








3/2 


9/2 


492.3817 




[0,1] 








3/2 


7/2 


475.5771 


[0.040656,0.999173] 


3 


1 


1/2 


7/2 


-915.7408 


[-0.999173,0.040656] 






3/2 


5/2 


413.6810 


[0.059441,0.998232] 






1/2 


5/2 


-872.0486 


[-0.998232,0.059441] 






3/2 


3/2 


336.9246 




[0,1] 





TABLE III: Hyperfine splitting (in MHz) and eigenstates for the ro-vibrational levels (v, L) with 
L = l,3 and v = 0, 1. All digits are converged. The relative theoretical accuracy on the frequency 



shifts, as well as on the smaller of the two coefficients 
few tens of kHz for the frequency shifts. 



ct,c£ 



, is 0(a 2 ). This corresponds to a 
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U 

91 = "1 


TT 

91=0 


91 = 1 


92 = -1 


( 2 ) X 

a g = o -2 

H HI ^ 

4 0) = o 


y 2 hi - 1 
af = 


„(2) _ V6r 

Jo) _ VI 
"0-3 


7T 

92 = 


4 0) = o 


tt Q — \/ ^a,0 
"0-3 


„(2) _ v^A 

4 0) = 


92 = 1 


„(2) _ v^A 
°g - — d q,0 

JP) _ ^3 
"O—3 


4 2) = f s q ,i 
4 0) = 


(2) X 
Oq = 0g i2 

4 0) = 



TABLE IV: Values of the coefficients a q for all combinations of the standard polarizations. 



L 


(ol||q(°)||ll) 


(ol\\qW\\il) 





0.7255 





1 


1.261 


0.7753 


2 


1.640 


0.8541 


3 


1.962 


0.9903 



TABLE V: Reduced matrix elements of the operators and for the transitions (v = 0,L) — > 
(V = l, L) with < L < 3, in atomic units. 



L 


A/ 


J 


J' 


7T7T 




<J + (T 





0.0000 


1/2 


1/2 


0.1754 


0.0000 


0.1754 




-50.7600 


5/2 


3/2 


0.0039 


0.0058 


0.0010 


2 


-1.2955 


5/2 


5/2 


0.1949 


0.0233 


0.1832 




1.9432 


3/2 


3/2 


0.1929 


0.0204 


0.1827 




51.4077 


3/2 


5/2 


0.0058 


0.0088 


0.0015 



TABLE VI: Average two-photon matrix elements S Q qi ,q 2 given by Eq. (|30p between the ro- 
vibrational levels (v = 0, L) and (v = 1,L) with L = 0, 2, in atomic units. A/ is the hyperfine 
shift of the transition frequency in MHz. 
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A/ 


(F: J) 


(F', j') 


7Y7V 




-|- 






—693.869 


(3/2,3/2) 


/ -\ 1 r\ o ln\ 

(1/2,3/2) 


1.028e- 


05 


1.534c 


-05 


2.607c- 


-06 


—689.945 


(3/2,5/2) 


(1/2,3/2) 


2.550e- 


06 


3.824e 


-06 


6.374e- 


07 


—684.601 


(3/2,3/2) 


(1/2,1/2) 


1.003e- 


05 


1.505e 


-05 


2.509c- 


06 


—680.678 


(3/2,5/2) 


/I 1 t-i 1 /*"i\ 

(1/2,1/2) 


1.118c- 


05 


1.677c 


-05 


2.794c- 


-06 


—645.591 


(3/2,1/2) 


(1/2,3/2) 


5.090e- 


05 


7.635e 


-05 


1.273e- 


05 




13/2,1/2) 


(1/2,1/2) 


4.229e- 


07 


0.000c 


-01 


4.229c- 


-07 


— 51.979 


(3/2,3/2) 


(3/2,1/2) 


1.329c- 


03 


1.993e 


-03 


3.322e- 


04 


— 48.056 


(3/2,5/2) 


(3/2,1/2) 


8.003e- 


03 


1.201e 


-02 


2.001e- 


03 


-10.348 


(3/2,3/2) 


(3/2,5/2) 


1.683c- 


02 


2.524e 


-02 


4.207e- 


03 


— 6.714 


(3/2,3/2) 


(3/2,3/2) 


1.853c- 


01 


1.281e 


-02 


1 . 789e - 


01 


—6.424 


(3/2,5/2) 


(3/2,5/2) 


1.842c- 


01 


1.122e 


-02 


1.786e- 


01 


— 3.702 


(3/2,1/2) 


(3/2,1/2) 


1.767c- 


01 


O.OOOe 


-01 


1.767e- 


01 


— 2.791 


(3/2,5/2) 


(3/2,3/2) 


1.122c- 


02 


1.683e 


-02 


2.804e- 


03 


2.487 


(1/2,1/2) 


(1/2,3/2) 


2.666e- 


02 


3.999e 


-02 


6.665c- 


03 


11.754 


(1/2,1/2) 


f 1 /<"i 1 1 n\ 

(1/2,1/2) 


1.767c- 


01 


O.OOOe 


-01 


1.767c- 


01 


12.325 


/I In O / O \ 

(1/2,3/2) 


(1/2,3/2) 


1.901e- 


01 


2.002c 


-02 


1.801c- 


01 


21.592 


(1/2,3/2) 


(1/2,1/2) 


1.333e- 


02 


2.000c 


-02 


3.333e- 


03 


37.929 


(3/2,1/2) 


(3/2,5/2) 


2.401e- 


02 


3.601e 


-02 


6.002c- 


03 


41.563 


(3/2,1/2) 


(3/2,3/2) 


2.657c- 


03 


3.986c 


-03 


6.643c- 


-04 


644.376 


(1/2,1/2) 


(3/2,1/2) 


4.229c- 


07 


O.OOOe 


-01 


4.229c- 


07 


654.214 


(1/2,3/2) 


(3/2,1/2) 


2.387c- 


05 


3.581c 


-05 


5.968e- 


06 


686.008 


(1/2,1/2) 


(3/2,5/2) 


3.637c- 


05 


5.455c 


-05 


9.092c- 


06 


689.641 


(1/2,1/2) 


(3/2,3/2) 


2.000c- 


05 


3.000c 


-05 


4.999e- 


06 


695.845 


(1/2,3/2) 


(3/2,5/2) 


4.102c- 


06 


6.153c 


-06 


1.026c- 


06 


699.479 


(1/2,3/2) 


(3/2,3/2) 


1.020c- 


05 


1.522e 


-05 


2.588c- 


06 



TABLE VII: Same as Table EH with L = 1. 
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A/ 


J) 


(F' : j') 


7T7T 


4- 


-\- 


+ - 




— 711.499 


(3/2,9/2) 


(1/2,7/2) 


6.739c- 


06 


1.011c 


-05 


1.685e- 


-06 


— 702.633 


(3/2,7/2) 


(1/2,7/2) 


4.143e- 


06 


5.629e 


-06 


1.329e- 


06 


—689.653 


(3/2,9/2) 


1 1 /ft c / <-* \ 
(1/2,5/2) 


1.179c- 


06 


1.768e 


-06 


2. 947c- 


07 


—680.787 


(3/2,7/2) 


(1/2,5/2) 


6.647e- 


06 


9.971e 


-06 


1.662c- 


-06 


—669.687 


(3/2,5/2) 


(1/2,7/2) 


1.073e- 


07 


1.610e 


-07 


2.684c- 


08 


G. A *7 O A 1 


to t n c / o\ 

(3/2,5/2) 


(1/2,5/2) 


1.457c- 


05 


2.030c 


-05 


4.417c- 


06 


—628.647 


(3/2,3/2) 


(1/2,7/2) 


1.764c- 


06 


2.647c 


-06 


4.411e- 


07 


— oOu.oOl 


to t n o /o\ 

(3/2,3/2) 


(1/2,5/2) 


3.055e- 


05 


4.582c 


-05 


7.637c- 


06 


-85.166 


(3/2,9/2) 


(3/2,3/2) 


O.OOOe- 


01 


0.000c 


-01 


O.OOOe- 


01 


— 76.301 


(3/2,7/2) 


(3/2,3/2) 


5.328e- 


04 


7.992e 


-04 


1.332e- 


04 


— 46.788 


(3/2,9/2) 


(3/2,5/2) 


3.324e- 


04 


4.986c 


-04 


8.310e- 


-05 


— 43.355 


(3/2,5/2) 


(3/2,3/2) 


5.742c- 


03 


8.613e 


-03 


1.436e- 


03 


— 37.922 


(3/2,7/2) 


(3/2,5/2) 


5.624e- 


03 


8.437e 


-03 


1.406e- 


03 


— 15.840 


(3/2,9/2) 


(3/2,7/2) 


4.070c- 


03 


6.106e 


-03 


1.018e- 


03 


-10.540 


(1/2,5/2) 


(1/2,7/2) 


2.677c- 


03 


4.015e 


-03 


6.691e- 


04 


— 7.438 


(3/2,9/2) 


(3/2,9/2) 


1.975e- 


01 


2.140e 


-02 


1.868e- 


01 


—6.974 


(3/2,7/2) 


(3/2,7/2) 


1.906e- 


01 


1.1 14e 


-02 


1.851e- 


01 


—4.977 


(3/2,5/2) 


(3/2,5/2) 


1.881c- 


01 


7.309e 


-03 


1.844e- 


01 


— 2.315 


(3/2,3/2) 


(3/2,3/2) 


1.922c- 


01 


1.345e 


-02 


1.855e- 


01 


1.428 


(3/2,7/2) 


(3/2,9/2) 


5.087e- 


03 


7.631e 


-03 


1.272e- 


-03 


11.306 


(1/2,5/2) 


(1/2,5/2) 


1.992c- 


01 


2.394c 


-02 


1.872e- 


01 


12.681 


(1/2,7/2) 


(1/2,7/2) 


1.999e- 


01 


2.499e 


-02 


1.874c- 


01 


25.971 


(3/2,5/2) 


(3/2,7/2) 


7.498e- 


03 


1.125e 


-02 


1.875e- 


03 


34.374 


(3/2,5/2) 


(3/2,9/2) 


5.538e- 


04 


8.308c 


-04 


1.385e- 


-04 


34.527 


(1/2,7/2) 


1 1 / n c / <-* \ 
(1/2,5/2) 


2.008e- 


03 


3.012c 


-03 


5.019c- 


04 


36.063 


(3/2,3/2) 


(3/2,5/2) 


8.616e- 


03 


1.292c 


-02 


2. 154c- 


-03 


67.012 


(3/2,3/2) 


(3/2,7/2) 


1.066e- 


03 


1.599e 


-03 


2.664c- 


-04 


75.414 


(3/2,3/2) 


(3/2,9/2) 


0.000c- 


01 


0.000c 


-01 


0.000c- 


01 


615.793 


(1/2,5/2) 


(3/2,3/2) 


2.202e- 


05 


3.304e 


-05 


5.506e- 


-06 


639.014 


(1/2,7/2) 


(3/2,3/2) 


9.467c- 


07 


1.420e 


-06 


2.367c- 


07 


654.171 


(1/2,5/2) 


(3/2,5/2) 


1.136c- 


05 


1.548e 


-05 


3.614c- 


06 


677.392 


(1/2,7/2) 


(3/2,5/2) 


2.496c- 


07 


3.745e 


-07 


6.241c- 


08 


685.119 


(1/2,5/2) 


(3/2,7/2) 


1.062c- 


05 


1.593e 


-05 


2.654c- 


06 


693.522 


(1/2,5/2) 


(3/2,9/2) 


2.124c- 


06 


3.186c 


-06 


5.310c- 


07 


708.340 


(1/2,7/2) 


(3/2,7/2) 


3.025c- 


06 


3.951e 


-06 


1.049c- 


06 


716.743 


(1/2,7/2) 


(3/2,9/2) 


9.039c- 


-06 


1.356e 


-05 


2.260e- 


06 



TABLE VIII: Same as Table EH with L = 3. 



L 


u 2ph (MHz) 


^2 P h (A*m) 





32 844 161.844 


9.128 


1 


32 798 213.622 


9.141 


2 


32 706 607.796 


9.166 


3 


32 569 919.581 


9.205 



TABLE IX: Spin-independent frequency and wavelength of the (t> = 0, L) — > (v = l,L) transitions, 



with < L < 3. They were calculated using the data of Refs. 
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